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We used a scries of CO^BOLD hydrodynamical model 
atmospheres covering stellar objects from vifhite dwarfs to 
red giants to derive theoretical estimates of the photo- 
metric and photocentric stellar variability in wavelength- 
integrated light across the Hertzsprung-Russell diagram. 
We validated our models against solar measurements from 
the SOHO/ViRGO instrument. Within our set of models 
we find a systematic increase of the photometric as well 
as photocentric variability — which turn out to be closely 
connected — with decreasing surface gravity. The esti- 
mated absolute levels of the photocentric variability do 
not affect astrometric observations on a precision level ex- 
pected to be achieved by the GAIA mission — with the 
exception of close-by giants. The case of supergiants re- 
mains to be investigated. In view of the ongoing debate 
about the photometric non-detection of p-modcs in Pro- 
cyon by the Canadian Most satellite we remark that we 
obtain a factor of w 3 in amplitude between the gran- 
ular background "noise" in the Sun and Procyon. This 
statement refers to a particular representation of tempo- 
ral power spectra as discussed in Sect. El 
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1. Introduction 

The presence of a small-scale time-dependent granulation 
pattern on the surfaces of late-type stars (see Fig. ^) leads 
to low-level temporal fluctuations in a star's total lumi- 
nosity, apparent position given by its photocenter, (see 
Fig. [SJ, and spectroscopically measured radial velocity. 
During recent years this kind of micro- variability got into 
the focus of research since it forms an important noise 
source in searches for solar-like oscillations and exoplanets. 
This is particularly the case for already operational or up- 
coming satellite mission like Most, Corot, and Kepler; 
granulation induced micro- variability might also become a 
limiting factor for precision astrometry missions like SIM 
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Figure 1. Typical granular intensity pattern from the red 
giant model sequence. Note the spatial scales. 



estimates of the convection-related photometric and pho- 
tocentric variability across the Hertzsprung-Russell dia- 
gram. To our knowledge the work of Trampedach et al. (1998) 
is the only previous example of such a theoretical effort 
who studied the brightness and radial velocity variability 
in the Sun, a Cen A, and Procyon A. 



2. CO^BOLD RAD. -HYDRODYNAMICS SIMULATIONS 

We used the radiation-hydrodynamics code CO'^BOLD 
(for further information about the code and applications 
see |Freytag et al. 2002| and |Wedemeyer et al. 2004| ) to con- 
struct a series of 3D Cartesian "local-box" model atmo- 
spheres. Table ^ summarises the model properties. We 
calculated three solar models which differ in numerical 
details. They are used as basic reference and for validat- 
ing our models by comparison with SOHO/ViRGO observa- 
tional data. The two models "Procyon A" and "^ Hydrae" 
have parameters close to the actual parameters of the stars 
they are named after. The other models are not intended 
to represent particular stars but are generic spanning a 



and GAIA (|Hatzes 2002IIGreen et al. 2003l[Sigrain et al. 20]Mfe ;e range primarily in surface gravity. Two models are 
IMatthews et al. 2004|l . We performed radiation-hvdrodvnamin^ended to investigate effects of metallicity. Note, that 
simulations for stellar atmospheres to derive theoretical the variation in effective temperature among the models 
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Figure 2. A statistical realisation of the motion of the 
photocenter in the plane of the sky based on data from 
the red giant simulation. The time between two diamonds 
amounts to 10^ s. 



is not very large, i.e. the set of models does not fully cover 
the parameter space. While not apparent from the table 
significant effort was invested to follow the evolution of 
the models over long periods of time to ensure statisti- 
cally representative results. 



3. From simulation boxes to stars 

The CO^BOLD simulations provide a statistically repre- 
sentative, rectangular patch (or "tile" ) of the emergent ra- 
diation field and its temporal evolution. In order to derive 
disk-integrated, observable quantities we have to extrap- 
olate this information to the whole visible stellar hemi- 
sphere. To this end we envision the stellar surface being 
tiled by a possibly large number of simulation patches. 
Being interested in integrated quantities we can ignore 
the spatial structure within in a patch and only need 
to consider averages of the emergent intensity Im over 
the patch surface as a function of time and inclination 
Am = cos(i?m). Due to limitations of computing resources 
the discretisation of the radiation field in solid angle is 
rather coarse, and usually we have only a total number of 
inclinations M (see Tab. of two or three available. 

Figure 13 shows an example of a resulting time series 
with M = 3. Clearly visible is the effect of limb-darkening 
and residual intensity fluctuations after averaging over the 
patch's surface. The crucial assumption now is that the 
size of a patch is large enough so that its emission can be 
considered as statistically independent of all other patches 
tiling the surface. Moreover, we assume that all patches 
share the same statistics given by the statistics of the sim- 
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Figure 3. Intensity time series for the three inclinations 
of model Sunl showing residual fluctuations and overall 
limb- darkening. 



rived expressions for the statistics of disk-integrated ob- 
servables which we summarise below. 



Statistics of the photometric variability. We characterise 
the photometric variability by the temporal power spec- 
trum of the relative fluctuations of the observable flux /. 
The expectation value of a frequency component denoted 
by f{i') is then given by 
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ulated patch. With these assumptions Ludwig (2004) de- 



Angular brackets denote expectation values, an asterisk 
the conjugate complex. The sums are discrete analogs of 
integrals over one half of the total solid angle where Wm 
is the integration weight. The numerator expresses that 
a frequency component of the resulting power spectrum 
is a Am weighted integral of the frequency components of 
the individual intensity power spectra. The denominator 
is essentially the square of the average observable flux. 
The spectral power density scales inversely proportional 
to the number of patches N tiling the visible hemisphere 
which is given by 

NA^2nR'^, (2) 

where A is the patch surface area and R is the stellar 
radius. R is assumed for the star the atmosphere model 
is associated with. In this paper we used plausible but 
to some extend arbitrary values for the stellar radii. For 
possible later adjustments of the radius one should keep in 
mind that the power of the photometric fluctuations scales 
as R~^ . The power spectra presented in this paper are 
given as spectral density normalised so that the integral 
between zero and the Nyquist frequency corresponds to 
the variance of the original signal. 

Statistics of the photocentric displacement. We charac- 
terise the absolute photocentric displacement along an ar- 



Table 1. CO^BOLD radiation-hydrodynamics model atmospheres: "Model" is a model's nickname used in this paper, 
Toff the effective temperature, logp the gravitational acceleration, R an assumed stellar radius (not intrinsic to the 
simulation proper), I the linear horizontal size of the computational box, the sound crossing time over i?™'^, i?™''^ 
the pressure scale height at Rosseland optical depth unity, (5/rms the relative spatial white light intensity contrast, Nobm 
the number of equivalent frequency points considered in the solution of the radiative transfer equation (RTE), M the 
number of inclinations used in the discretisation of the RTE, "Modelcode" an internal identifier of the model sequence. 
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bitrary axis a;ph by its standard deviation (Jxph given to 
good approximation by 
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where I is the Hnear patch size, ct/ the standard deviation 
of the observable flux, and (T/„^ the standard deviation 
of the intensity at inclination The first equality in 
Eq. (O expresses a close relationship between the level of 
fluctuations in the photocentric position and the bright- 
ness. This relation carries also over to the temporal power 
spectrum of Xpii{t) which is up to a factor R^/6 identi- 
cal to the power spectrum of the brightness fluctuations. 
Remarkably, the last equality in Eq. lO expresses that 
Cxph is not dependent on stellar radius since only quan- 
tities intrinsic to the simulation enter. The flrst equality 
might suggest the contrary. However, cr//(/) scales in- 
versely proportional to the radius so that the radius factor 
cancels. The similarity in the power spectra has the con- 
sequence that the characteristic time scale for the fluctu- 
ations in photocentric displacement and brightness is the 
same, namely the evolutionary time scale of individual 
convective cells. The similarity in the power spectra does 
not imply that the fluctuations are correlated; it can be 
shown that photocentric and brightness fluctuations are 
perfectly uncorrelated. 

4. Comparison with VIRGO observations 

FigureElshows a comparison between the disk-integrated, 
photometric fluctuations derived from our three solar mod- 
els, and observational data from the ViRGO instrument 



on board the SOHO satellite. Our focus is the high fre- 
quency region of the solar signal in which granular contri- 
butions dominate. The ViRGO power spectrum has been 
calculated from (level 2) time series data provided^ by 
the Virgo team. Actually chosen was one year of data 
(1996.5-1997.5) close to solar minimum activity to min- 
imise the possible contribution of activity related vari- 
ability. Moreover, the plotted power spectrum is based on 
data of the green channel of the ViRGO three-channel sun- 
photometer (SPM, see lFrohlich et al. 1997|l . It has been 
converted to white light fluctuations by matching (by shift- 
ing in power) a corresponding power spectrum based on 
Virgo PM06V-A absolute radiometer data in the fre- 
quency range 0.3 to 2.0 mHz. Instead of using the PM06V- 
A power spectrum directly this rather involved proce- 
dure was necessary since PM06V-A and SPM power spec- 
tra deviate substantially in the high frequency region. 
The authors found little information about this mismatch 
in the literature. However, it appears to be consensus 
that the SPM spectrum reflects the actual solar behaviour 
(|Frohlich et al. 19971 lAndersen et al. IQQ'Sji . in particular 
showing the steep (roughly as i/""^) decline at the highest 
frequencies. 

We find a reasonable agreement of the continuous back- 
ground signal between model predictions and observations 
between 0.04 and 8.0 mHz. In particular, the background 
in the p-mode frequency range is matched quite well. The 
differences among the models can be taken as an estimate 
of the involved numerical and statistical theoretical uncer- 
tainties. We emphasise that the also visible p-mode peaks 
of the models are not expected to match directly the ob- 
servations since the resonance and excitation properties 
differ widely between the box models and the Sun. 



^ |http://www.pmodwrc.ch/pmod.php?topic=tsi/virgo/proj-space-virgo#'\ 
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Figure 4- Comparison of power spectra of disk-integrated 
photometric fluctuations between three solar simulations 
and observational solar data from SOHO/VIRGO. Note 
the steep decline in power in the range of the p-mode fre- 
quencies. 



5. Trends in the photometric variability 

Figure El shows a comparison of the brightness fluctua- 
tions for a sequence of models spanning a range between 
DA white dwarfs with convective outer envelopes and red 
giants. The power spectral density Pi, is plotted as vPi, 
which makes it independent of the unit in which the fre- 
quency is measured. This facilitates the intercomparison 
among the models where we scaled the frequency with 
the sound crossing time over a pressure scale height at 
the surface (tc). As evident from Fig. [S] all spectra show 
a similar shape, and are essentially located in the same 
scaled frequency range. We find a systematic increase of 
the photometric variability towards giants. Note, that in 
fact the square root of the power is plotted in Fig. [S] I.e. 
our models predict an increase in the amplitude by a factor 
of « 1000 between the white dwarf and red giant model. 

Our result is in marked contrast to the modelling of 
Trampedach et al. (1998) who did not find an increase in 
the granular photometric signal comparing the Sun and 
Procyon A. However, they pointed out that their time se- 
ries might have been not long enough to provide sufficient 
statistics and coverage of lower frequencies. In view of 
the ongoing debate about the photometric non-detection 
of p-modes in Procyon by the Canadian Most satellite 
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Figure 5. Spectral power density of the brightness fluctua- 
tions as a function of scaled frequency v -Tc- Diamonds de- 
limit frequency ranges in which p-modes have been detected 
in radial velocity (Procyon data were taken from Brown et 
al. 1991, ^ Hydrae data from Frandsen et al. 2002); ob- 
served p-mode frequencies are located in regions with large 
background power in the models. Note the systematic in- 
crease in power with decreasing gravity. The behaviour at 
scaled frequencies greater than 300 is dominated by the 
numerics and should be discarded. 
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(|Matthews et al. 2004l|Christe nsen-Dalsgaard & Kjeldsen 
we stress that we obtain a factor of « 3 in amplitude 
between the background signal in the Sun and Procyon. 
This refers to a comparison in a representation like Fig. [3 
where differences in power can be essentially described by 
uniform vertical shift. In a representation in absolute fre- 
quencies the power ratio would be frequency-dependent. 

By comparing models of similar effective temperature 
and surface gravity, but different metallicity we find a 



Figure 6. Spectral power density of the brightness fluctua- 
tions for two models around Toff = 5000 K of solar (filled) 
and 1/100 solar metallicity (dotted). To facilitate compar- 
ison the X-axis includes the same range of frequencies as 
the one in Fig. \^ 

2004| ) 



decrease of the brightness fluctuations with decreasing 
metallicity, see Fig. El This is in line with expectation that 
the higher densities encountered at optical depth unity in 
metal poor models (due to lower overall opacity) leads to 
smaller convective fluctuations and consequently smaller 
brightness fluctuations. 



Figure 7. Photocentric displacement as a function of sur- 
face gravity. The diamond indicates the metal poor model. 



6. Trends in the photocentric displacement 

We find an almost perfectly linear relationship between 
the standard deviation of the photometric displacement 
and the surface gravity of the models (see Fig. Ej). One 
has to keep in mind, however, that there will be some 
scatter around this line when considering atmospheres of 
markedly different effective temperature or chemical com- 
position from the ones considered here. The diamond de- 
picting the metal poor model already indicates this. 

Instead of the basic atmospheric parameters (Toff, log g, 
chemical composition), two parameters more closely re- 
lated to the convection pattern are perhaps physically bet- 
ter suited to describe the functional behaviour of the pho- 
tocentric displacement — namely the relative spatial in- 
tensity contrast (J/mis and the typical size of a convcctive 
cell. The size of granular cells is related to the pressure 
scale height at optical depth unity i?^"''^( |EY'eytag 20011 ). 
Figure|H|depicts the outcome of a test of a relation between 
photocentric displacement and the product SlmisHp^'^^' in- 
deed, all models including the metal poor one now follow 
the same trend. This also indicates that the proportional- 
ity to g~^ seen in Fig.Qis primarily reflecting the system- 
atic increase of granular cell size with decreasing gravita- 
tional acceleration. 

For our most extreme giant model with logg = 1 we 
find a standard deviation of the photocentric position of 
3 • 10"''' AU or 0.3 mas/Z? [pc] where D denotes the dis- 
tance. On a precision level expected to be reached by 
the GAIA astrometry such a variability does only affect 
the achievable precision for close-by giants. The situation 
for supergiants remains to be investigated; a linear ex- 
trapolation of the curves shown in Figs. |7| and |H| towards 
lower surface gravity is unlikely to give a reliable esti- 
mate of their photocentric variability since sphericity ef- 
fects render convection in supergiants markedly different 
( [Freytag et al. 2002| ) from convection in Cartesian geom- 
etry as presented here. 



Figure 8. The displacement as function of relative spatial 
intensity contrast times surface pressure scale height. 



7. Concluding remarks 

We found systematic changes of the convection-induced 
fluctuations in brightness and photocentric displacement 
with stellar parameters. The result is a nice example for an 
application of hydrodynamical model atmospheres (in this 
respect see Ludwig & Kucinskas, this volume). We would 
like to emphasise again that we have been considering the 
contribution of granulation to the stellar variability only. 
At highest temporal frequencies it is the dominant con- 
tributor to the variability. However, at lower frequencies 
variability related to magnetic activity, spottedness, and 
rotation dominates, and the level of the fluctuations is 
likely to be significantly larger than the convective contri- 
bution. 
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